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Abstract 

We consider induced pair production in an external field at finite temperature. 
One-loop correction to the Green function of a meson is calculated semiclassically 
within the framework of saddle-point analysis of Schwinger proper time integrals. 
This correction appears to be exponentially small in terms of inverse temperature 
dependence. Low-temperature limit is shown to be in full agreement with previ- 
ously obtained zero-temperature results. The corrections in the low-temperature 
limits are estimated up to the leading exponential and pre-exponential terms. 
Comparison is made to earlier calculations of vacuum decay. 
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1 Motivation 



Spontaneous processes of particle production in field theory [JJ (also known as 
Schwinger processes) or string/brane production [2] in string theory in external 
fields have long been studied. Production of e + e~ pairs by a constant electric field 
is the archetypal example for the wide class of these non-perturbative phenomena. 
They can generally be characterized by the essentially non-analytic behaviour of 
observables in the external field in the weak field limit, that is, by the presence of 
e~ "B-like terms. A closely related class of phenomena is known as vacuum decay 
processes OH]- One of the similarities between Schwinger pair production and 
vacuum decay is that they both can be described in terms of the semiclassical 
approximation to the tunneling problem in quantum mechanics. That is, the 
leading probability or another observable is usually organized as e~ SE , where Se 
is some action on some classical configuration. 

The process being spontaneous means the initial state must be a vacuum 
state. A generalization of the Schwinger phenomena to the processes in which a 
non-zero excitation is contained in the initial state (i.e. some particle is present) 
is referred to as induced Schwinger process. Induced processes of vacuum decay 
have also been known for quite a long time [5j Ej. 

Recent progress in understanding induced brane production [2] at zero tem- 
perature has lead us to the following question: how would finite temperature 
influence the dynamics of brane production/decay? This would be of great im- 
portance for cosmology. Brane induced decay is apt to be viewed upon as false 
vacuum-decay in a higher-dimensional theory. This task, when simplified down 
to field theory level, may be presented as particle decay in an external field. Re- 
cently decay of a magnetic monopole was calculated by the authors of the present 
paper within this string-motivated paradigm, the same was done for a Thirring 
model meson decay in 2D theory at zero temperature [7J. Here the results of [7J 
are generalized towards the case of finite temperature. 

Studies of spontaneous Schwinger pair production at finite temperatures have 
a long history. A number of papers have been produced during the last three 
decades on the spontaneous process of pair production in one-loop approxima- 
tion. Not claiming to have made a full review in the least part, we cite just a 
few of them \§\ [TU] . A modern picture of one- loop thermal results is reflected 
in [UJ. Two- loop results are available as well [12]. On the side of vacuum de- 
cay, the finite-temperature effects have been thoroughly studied by [TS] for the 
spontaneous case, yielding both the exponential and the pre-exponential contri- 
butions to the decay rate. On the side of the induced decays, the pre-exponential 
factor was first calculated in |14j . 

However, no results dealing with temperature corrections to induced decays 
are available so far. This was one of the main motivations for writing this paper. 
While figuring out the simplest one-loop correction to the propagator of a scalar 
particle due to temperature and external field, semiclassical treatment is applied 
which is very close to the world-line techniques by Dunne et al. |15|. I16j . The 
semiclassical approach to Schwinger processes has been suggested since Popov's 
papers, see e.g. [TT] . 

This article is organized as follows. In Section 12.11 a brief reminder of Green 
function techniques is given in the finite-temperature field theory, and the cor- 
rection to the meson Green function is calculated, giving a universally valid 
expression (in terms of any regime in = y). In Section [2.21 its asymptotics are 
studied for (3 — * oo, and comments are made on the opposite limit in Section 
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12.31 In Section 12.41 the problem of leading and sub-leading asymptotics in semi- 
classical calculations, as well as on a relationship between loop resummation and 
multi-instanton resummation is briefly discussed. We conclude in Section [3l 



2 Schwinger Processes at Finite Tempera- 
ture 

2.1 General Techniques 

A simple cubic interaction of a charged scalar eft and a neutral \ scalar in a 
two-dimensional theory is considered, 



C = \\D^\ 2 - ^/i 2 H 2 + \{d^ X f - \m 2 X 2 + A#* X , (1) 

where the covariant derivative given as = d^ + ieA^. The masses of the fields 
being \i and m are first kept arbitrary, but after proceding to the semiclassical 
approximation it will be assumed that - < 1. This situation is known in 

A* r—n 

vacuum decay terms as "an almost spherical bubble" and is used, e.g. in [14j . 
The charged field interacts with an external Abelian field An- 

One can show [18] that in the coordinate representation a free Green function 
at zero temperature for a particle of field x with zero charge is 

1 da. ■ 2 i{y- x ) 2 

where we have omitted the pole prescription is. For the charged particle <f> in 
the constant external field An = (0,Exq) it becomes 

G, \ 1 f ddo, im 2 a — (y— x) 2 coth(ea)— (Vl~ Xi)(yo+Xn) /o\ 

<t>( x >y = 7T~\2 / ^T7 — V e 4 4 > (3 

(47r)" ! ,/ smh(ea) 

2 

where e = e£7; the field is considered to be far below the Schwinger limit S> 
1. We shall refer to representation (|3|) as Schwinger parametrization, and the 
variable a — Schwinger parameter. When the temperature is not equal to zero, 
the Green function is periodic in the Euclidean time with a period /3 = ^, 
hence its generalization is organized as 

GJxil) = / £ ^ a c im 2 a — 1 ^\(yo-x +nl3) 2 + (y 1 -x 1 ) 2 ] coth(ea) ^ 

Jo (4vr) 2 sinh(ea) 
x e - j(m-^i)(2/o+a;o+n/3) 

(4) 

The sum over n naturally appears since one should take into account all the 
equivalent positions separated by (3 in Euclidean time, as explained in the cited 
paper. 

Let us consider the one-loop perturbative correction to the Green function of 
the uncharged particle due to the cubic interaction mentioned above 

6G X (0, z) = A 2 J d 2 xd 2 yG x (0, x)G 4> (x, y)G (j) (y, x)G x (y, z). (5) 

We think of a one-loop diagram, represented in configuration space in Figure 
[TJ We stress here that our adherence to configuration space representation of 
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Figure 1: One- loop vacuum polarization correction to the propagator of field \ i n an 
external field at finite temperature. Geometrical meaning of the saddle-point value a 
for the Schwinger parameter a is illustrated; 0, x, y, z correspond to the similar variables 
in ([HD; </>, X denote propagator type. 



vacuum polarization is not an incidental or technical detail of calculation, but 
is rather of conceptual meaning. Namely, as we have shown in our previous 
paper [7], the saddle-point values of on correspond directly to the geometric 
parameters of the classical sub-barrier (Euclidean) trajectory in configuration 
space. The other reason to keep the position space is the direct relationship of 
the spatial configuration shown in Figure[T]to a vacuum bubble with external lines 
attached to it [H] describing induced vacuum decay on a compactified manifold. 
The correction to Green function of the neutral \ field becomes then 



A 2 e 2 \ f°° da± dcti da.3 da^ d 2 x d 2 y 



sinh(eoi) sinh(eo2) 03 «4 

p,q,k,n 

i(x +kf3) 2 _ i(z -y Q +n(3) 2 _ _ i{zi_-y{) 2 ijJ3_ _ . . 
x e 4a 3 4a 4 4a 3 4a 4 ~r 2 W 1 X1 ' l - P q > 

x e _ T^ yo ~ xo+p ^ 2 + ( yi ~ X1 ^ coth ( ea l) — ^[{yo-xo+ql3) 2 + {yi-x 1 ) 2 ]coth(ea 2 ) 

(6) 

Note that due to temperature, Lorentz-invariance is explicitly broken in this 
expression. We may relate coordinate representation of the Green function cor- 
rection to 

5G(n,ki) = J ' SG(0;z , Zl )e- iuJnZo - iklZl dz dz x , (7) 

where oj n = ^ J w L - Calculation of this correction corresponds, as usually, to the 
shift of the Green's function pole 

— * % 

G(n ' k) = cj^ + k 2 1 +m^ + AP(n,k i y (8) 

however, now the shift is not Lorentz-invariant but rather depends on both n and 
k separately. Whatever complicated expression for the variation of the Green 
function one would obtain, physically relevant information is contained in the 
pole shift and rescaling of the wave function. Further on, only the pole shift will 
be considered. The mass shift is the value of the M 2 calculated at the point 
corresponding to the pole. For the zero temperature theory that means that the 
shift is determined only by the value of bare mass (and a scale). Since in this 
case M 2 (k 2 ) depends only on k 2 = k\ — k 2 , one should take k 2 = m 2 . For a 
theory in a compactified Euclidean space-time it is quite obvious that the shift 
should depend on n as M 2 (n, k±), where n and k\ such that u; 2 + k\ + m 2 = 0. 
For the n-th Matsubara mode the pole shift is approximately related to a Green 
function variation as 
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5m r 



in 



+ k 2 + 



47r 2 n 2 



P 2 



SG(k u n) 



(9) 



i+fca- 



Moreover the on-shell condition does not have a solution for arbitrary k\ since 
uj n is determined by the discrete variable n. Note again that generically every 
single mode is renormalized in its own way. One can easily understand that due 
to Lorentz symmetry violation by compactifying the Euclidean time direction 
there is no invariant mass anymore. Rather, if some physical quantities related 
to thermal theory are of interest, the partition function Z(f3, fi, . . . ) has to be 
calculated, at some values of chemical potential \x and other external potentials. 
Then some statistically averaged reasonable quantities have to be founcQ, say, 
rij = J^-, which is the equilibrium concentration of particles of the i-th type. 

Therefore, is understood merely as a convenient way of writing down the 
propagator variation. However, there is a range of parameters within which it 
is still possible to preserve the meaning of this quantity as the mass shift of the 
particle. This range is the limit of small temperatures or large (3. In this case 
on-shell condition can be solved even for k\ = 0, since for sufficiently large f3 the 
value of ^ differs fr om an integer slightly. So one can treat the imaginary part 
of the mass shift for such n and k\ = as the decay rate of the particle in the 
external field with a temperature not equal to zero. 

Evaluating elementary integrals, one gets a formal expression for the mass 
shift 



5m 



x 
x 



X 2 



+oo+i0 



g e 3/2 



Jn+r+s 



o 



y / sinh(ai + 02) cosh(ai — 02) 

4ir 2 i (r tanh(ai)— s tanh(Q2)) 2 sinh(2ai ) sinh(2a2) 
g e/3 4sinh(cti+a2) cosh(ai— 012) y_ 



r 2 tanh(ai)+s 2 tanh(«2) + (01+0:2) + 



tf}' 2 coth(«i)+coth(a2) 



(10) 



By using the approximate on-shell condition n ~ & and Poisson resummation 
formula 



where 



f(k) 



f(t)e~ lkt dt, 



this expression can be written down in two equivalent representations. Now 
and further we retain its imaginary part solely, rather then the full (possibly 
divergent) pole shift in the propagator. These two representations are 



r = Im — — 



A 2 ^ f da X da 2 e — (Q1+Q2) " — coth(a 1 )+coth(a 2 ) + 0) 



l fl e 3//2 J ^Jsmk{a\ + a 2 ) cosh(ai - a 2 ) 

(11) 

and 



: The authors are extremely grateful to Professor H. M. Kleinert for a discussion on this point. 
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ifJ. 2 , s im 2 1 ie/3 2 s 2 . 

A 2 ^ f d ai da 2 e~ {ai+a2) ~ — cothCaO+coth^) " TT-^ 

Im — > / — — - , 

me J Sinn fti + ctz) 

s=—oc 

(12) 



where 



A 



so 



sinh(ai + 02) 
cosh(ai — 02) ' 

m(3 1 



2ir tanh(a2) coth(ai) + 1 

The sums above can formally be converted to Jacobi theta- functions, however, 
that would not be of much use, since the integrals over Schwinger parameters 
would then get out of feasibility. On the contrary, one can do the integrals in 
Schwinger parameters by saddle-point method for each term in the sum, provided 
saddle-point works at all. Then in principle one could try to do the sum exactly. 
When dealing in such way with 1-dimensional saddle-point integrals 



& ifn(z) 



it will be necessary to restrict the domain of applicability of this approximation 
by imposing the conditions in saddle-point values z = z n for the nth function /„ 

mi-. 



\fo(zo)\ 3/2 ^\f{'m, (13) 
l/i^o)| 3/2 » /f (zi), (14) 

Im[/ (io)] > MAM- (I 5 ) 

Condition (fT3|) ensures the possibility of doing saddle-point approximation for 
the zero mode. It checks whether the next-to-leading order terms in the decom- 
position of /o may be neglected. Condition (|14|) provides the same check for 
the first mode. To ensure dominance of the zero mode over the first one, we 
impose (fT5j) , 

For multi-dimensional expressions the criteria of saddle-point method appli- 
cability become more complicated. Namely, instead of (|13p or (|14p one must 
require that 

d 3 f 1 

" -PuP jm Pkn^=== « 1 (16) 



OZiOZ-jOZk V^^raAn 

where A,- arc eigenvalues of second derivatives matrix q z q z 

-, and Pij is diagonal- 

isation matrix for q z .q z . , summation implied over i,j,k,l,m,n. Further it will 
be examined whether these conditions are satisfied for a particular saddle-point 
function under consideration. 
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2.2 Limit (3 — > oo 



The expression (I12p seems to be an appropriate representation for the V in the 
case (3 — > oo, since even the naive condition of the saddle point method applica- 
bility fails for pip , namely the factor in the exponent becomes small. So, let 
us take (|12p and make sure that it indeed corresponds to the low-temperature 
limit. It is supposed that the saddle-point is a symmetric point a\ = «2 = ex. 
This is possible due to having particles of identical mass in the internal lines. 
Since the saddle-point values of Schwinger proper times, as shown in [7J, corre- 
spond to the geometric parameters of the classical Euclidean loop configuration, 
only a symmetric configuration is expected to be the physically relevant solution 
of saddle-point equations. 

The full "decay width" is a sum over Matsubara contributions 

oo 

r= J2 r «- 

n=— oo 

It is expected that the higher the mode, the more suppressed it is. It will be 
shown below by means of saddle-point integral that this is indeed true for the 
zeroth and first modes. One can trivially see that the zero Matsubara mode 
exponential is 

fi 2 m 2 tanh(o) 
Jo = — 2a 



e e 2 ' 

identical to that of [7|. This function is extremized for a = 6tQ given by 

cosh(a ) = — • 

For simplicity the case of a very light external particle ("almost spherical bub- 
ble") will be considered, that is ^ <C 1. This is indeed the case of interest, as for 
a particle with m > 2/i the process will become perturbatively allowed. Checking 
the conditions of saddle-point applicability (fl6l) . one gets 



< 1. 



This condition is satisfied in our setting, due to having a sufficiently small field 

2 

— 3> 1, imposed from the very beginning for the applicability of the saddle-point 
method (there always must be a significant exponential suppression). Thus the 
leading order (zero dual Matsubara mode) contribution to the sum (|12|) is 



JL :e-^). 



yMetjj didjfo 

In the leading order in both small parameters ^ and one obtains 

det didjfo = -^r, 
hi e 2 

hence 

A 2 V 

e" — , (17) 



m 2 fj, 
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in agreement with [7]. The next-to- leading term is given by the modes with 
s = ±1. In the limit e/3 2 3> l,/i>m one obtains 



/ ±1 - ■ fl 



e 4 sinh 2a 

That amounts to the saddle-point equations 

cosh 2a = -z-zo, 

solved by 

iri 

in the said approximation. Collecting all the terms, one gets the first correction 



r ±l~ e 4 - ( 18 ) 



A 2 _M! 
me 2 /3 2 

It can be seen that the dependence on a temperature is essentially non- 
perturbative. 

Upon calculating the second and the third derivatives of f±i, one gets the 
following inequality as the condition for saddle-point method applicability by 
evaluating ([To]) 

16u 2 , , 

<K h (19) 

which can be rewritten as 



.2X1 



p»1{t) ■ < 20 > 

Obviously this condition is fulfilled provided that the temperature is high enough. 
The value of the exponential f±\ on the saddle point in the leading asymptotics 
is 

/±i(«o) = ^/3 2 e- 

In the aforementioned limit the s-th mode will be suppressed like e~ e ^ s2 . The 
condition for the effective suppression roughly is 

e 

which is satisfied for (3 — * oo limit. Thus one makes sure that the chosen form of 
the series (fT2|) indeed suites low-temperature region description. 



2.3 Limit (3^01 

Naively, it seems that one can easily perform similar calculations for the oppo- 
site case of extremely high temperatures. Just using the dual representation for 
the decay rate (fTT|) and applying the saddle point method one could get the 
answer. However, this is not necessarily so for several reasons. In both rep- 
resentations for r the pole has already been chosen by setting k\ = 0. It is 
easy to see that it is impossible to set kf = and find such an integer value n 
which would satisfy the on-shell condition ^-^r — V = for small /3. So, in 
order to fulfill the on-shell condition and, thus, to find the proper expression for 
the rate in the limit of high temperature, one should keep k 2 non-zero, which 
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reflects violation of Lorentz symmetry by introducing a temperature. It is also 
worth mentioning that in this case, since one cannot get rid of the dependence 
on n, it is absolutely necessary to consider the renormalization for each Matsub- 
ara mode separately. The result certainly cannot be interpreted then as a high 
temperature decay rate of the initial particle. It is just some correction to the 
Green function, the physical meaning of which is not well defined. It may be 
meaningful in a compactified theory rather than in a thermal one. Intuitively 
one expects that at high temperatures tunneling term e~ So (instanton contribu- 
tion to the semiclassical expression) will be dominated over by the over-barrier 
term e~^ £: ("sphaleron" contribution). However, in the saddle-point analysis per- 
formed by us, the "sphaleron" term has not appeared in the limit 3 — ► 0. Thus no 
final judgement is passed on the applicability of saddle-point method at 3 — > 0, 
but it is doubtful that it can work as directly as it has worked for 3 — » oo. 



2.4 Resummation Hierarchy 

The values of a given above are, of course, not unique. The true solution to the 
saddle-point equation is a series of roots like, say, a = ±i arccos P + 2mn, n £ 
Z. Therefore, an additional resummation to include all these roots is to be per- 
formed in principle. It is not necessary from a practical point of view, the terms 

in the series being suppressed by the factor e , however, this resummation 
is noted in order to stress the similarity of the saddle-point configuration, on 
which the 1-loop integral is essentially calculated, and the world-line instanton 
configuration, proposed in the semiclassical approach by Dunne et al. |15} I16j. 



3 Conclusion 

An example of a calculation of the one-loop thermal corrections to the propagator 
of a neutral scalar particle interacting with a charged one in an external field has 
been given. We have found the LO thermal correction to decay width in the 
semiclassical limit for 3 — > oo, in the case of an "almost spherical bubble", i.e. 

2 

— <C 1, far below Schwinger limit — 3> 1, up to its preexponential factor 



A 2 
me 2 3 2 



5r ~Z^2 e 4 - (21) 



Thus the result of [7] has been generalized towards thermal/compactified back- 
grounds. Of peculiar interest would be extension of the presented semiclassical 
techniques to strings and branes in thermal backgrounds. Our result supports 
that of Garriga [13]. Namely, the technique of Garriga does not intercept any 
correction for vacuum decay in two dimensions at low temperatures, whereas we 
give an estimate of this correction, and point out that due to its rapid decrease 
it can't have been noticed in the framework of Garriga's method. 
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